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SUMMARY 


We consider the modeling and active damping of an elastic beam using distributed 
actuators and sensors. The piezoelectric ceramic material (PZT) is used to build the 
actuator. The sensor is made of the piezoelectric polymer polyvinylidcne fluoride 
(PVDF). These materials are glued on both sides of the beam. For the simple clamped 
beam, the closed loop controller has been shown to be able to extract energy lrom 
the beam. The shape of the actuator and its influence on the closed loop system 
performance are discussed. It is shown that it is possible to suppress the selected 
mode by choosing the appropriate actuator layout. It is also shown that by properly 
installing the sensor and determining the sensor shape we can further extract and 
manipulate the sensor signal for our control need. 

1 INTRODUCTION 

There has been an increasing interest in the control of large space structures 
and flexible structures in recent years. These structural systems are usually large in 
size, light in mass and hence weakly damped. In order to achieve vibration suppression 
and precision pointing, it is necessary to introduce artificial damping to such systems. 
One approach is passive damping by adding the minimum weight of damping material 
to the effective locations on the structure. Another way is to use external mechanisms 
with feedback of the systems’ state or output to counteract the undesired mot ion 
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modes. In modern structural engineering, active feedback controls to stabilize the 
structure are preferred. 

Proper modeling is essential to control system design and to avoid spillover due to 
the infinite dimensional nature of these systems. Here we consider the beam model 
as part of the structure and study its modeling and active damping. 

The actuator considered here is a distributed one made of piezoelectric ceramic 
material which is glued to the beam. Its constitutive property, i.e. its strain and stress 
relation, is influenced by the external voltage applied to it. Bonding or embedding 
segmented elements of this material in a structure would allow the application of 
the localized strain to be transferred to the structure whose deformation can be 
controlled. Under proper bonding conditions, the coupling between the actuator 
strain and the beam strain can be determined to implement the control mechanism. 
In [6] [2] [7] active vibration control is described using spatially distributed actuators. 
The PVDF sensor is bonded to the beam in a similar way. The output voltage 
is a functional of beam curvature. Unlike the conventional point sensor, this is a 
distributed one. Cudney [5] provides some detailed explanation of the nature of 
piezoelectric actuators. 

We first discuss the modeling of the beam and the distributed actuator. A static 
model of the actuator coupled into the structure is developed. We have developed 
the beam model from the Euler- Bernoulli model with rotational inertia added. Next 
the sensor model is addressed. We then discuss the controller design using Lyapunov 
methods. We finally investigate the actuator and sensor shapes and their impact on 
the system elastic modes. 


2 SYSTEM MODEL 

One approach to build the desired actuator is to take advantage of the special con- 
stitutive properties of certain materials. The actuation is due to the property change 
under certain stimulation other than the external actuation force. Such materials 
are the so called smart materials. Once properly embedded into the structure the 
induced actuation will produce bending or stretching or both to control the structure 
deformation. One of the advantages of using smart materials as actuators and sensors 
is that the structure will not change much. 

Piezoelectric actuators were used as elements of intelligent structures by Crawley 
and de Luis [4]. Bailey and Hubbard [1] have used PVDF actuators to control the 
vibration of a cantilever beam. The control voltage applied across the PVDF is the 
sign of the tip rotation velocity multiplied by a constant. 

Figure 1 shows the structure of the beam with both the sensor and the actua- 
tor layers glued together. In this figure, h stands for the thickness of the different 
layers of the beam. The subscripts s, b and a denote sensor, beam and actuator 
respectively. The constitutive law for piezoelectric materials has several equivalent 
forms. The stress -strain relationship for the piezoelectric material is similar to that 
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Figure 1: The composite beam 

of thermoelastic materials, with the thermal strain term replaced by the piezoelectiic 
strain A. The constitutive equation of the actuator is given by 

a = E a {e- A) (1) 

where A is the actuation strain due to the external electric field, and e is the strain 
without external electric field. E a is the Young’s modulus of the actuator, cr is the 
stress of the actuator. The actuation strain is given by 

A(a :,t) = ^-V( Xl t) (2) 

fl>a 

where d^x is the piezoelectric field and strain field constant. V(x,t) is the distributed 
voltage. The strain has two effects on the beam. One effect is that it induces a 
longitudinal strain t[ to insure a force equilibrium along the axial diiection. This 
steady state value of Si can be derived by solving a force equilibrium equation. The 
other effect is that the net force in each layer acts through the moment arm with the 
length from the midplane of the layer to the neutral plane of the beam. The resultant 
of the actions produces the bending moment. Taking a similar approach as in [1] the 
actuation moment can be expressed as 

M a = ATA(M) (3) 

where K a is a constant depending on the geometry and the materials of the beam. 
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We model the beam with linear bending, no shear but with the rotational inertia 
included. This is more accurate than the Euler-Bernoulli beam model. The Euler- 
Bernoulli beam model 


+ P A ef ‘ w }*’ ,) = 0 


dx 4 


dt 2 


is often used to describe the beam dynamics because of its simple form and ease for 
analysis. E stands for Young’s modulus and I stands for moment of inertia. This 
equation can be rewritten as 


where 


* | EI 9 -^!} + = 0 


dx 


„2 I 


dx 2 


dt 2 
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The bending moment of the composite beam without actuation is 


( 6 ) 


M b = El 


d 2 w(x, t) 

dx 2 


The Euler-Bernoulli model is a linear model without accounting the rotational 
inertia and the shear effect. It is easy to see that during vibration the beam elements 
perform not only a translational motion but also rotate. The variable angle of rotation 
which is equal to the slope of the deflection curve will be expressed by dw/dx and 
the corresponding angular velocity and angular acceleration will be given by 


d 2 w 

dxdt 


and 


d 3 w 

fcdt 2 ' 


( 8 ) 


Therefore the moment of the inertia forces of the element about the longitudinal axis 
will be 

d 3 iv 


-pi 

The equation with rotational inertia is [8] 


dxdt 2 


( 9 ) 


d 4 w(x,t) d 2 w(x, t) d 4 w(x,t) 

~~a^- + pA ~Mi — = 


( 10 ) 


where 

pA = p a A a + pbAb + p s A s .. (II) 

We take this equation as our beam model under consideration. It falls in between the 
Euler-Bernoulli beam and the Timoshenko beam. 

The total bending moment with actuation is 


M = M b + M a . 


( 12 ) 
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Substituting Equation (12) into (10), we have 


d 4 w(x,t) d 2 \(x,t) A d 2 w(x,t) d 4 w{x,t) 

EI ~d. — ~ Ka ~d^~ + pA ~^e~ " pI ^w~ ~ °- 

Considering the actuation strain and the applied voltage, we get 


El 


d 4 w(x,t) d 2 V(x,t) d 2 w(x,t) d 4 w(x,t) 


dx 4 

with boundary conditions 


— c 


dx 2 


+ pA 


dt 2 


-pl- 


dx 2 dt 2 


= 0 


El 


El 


in(0, t ) 

<9u>(0, t ) 

dx 

d 2 w(L , t) 

fa 2 

d 3 w(L, t ) 

dx 3 


= c. 


0 

0 

c.V{L, t) 
dV(L, t) 


dx 


where 


i J 1 a • 

ha 


( 13 ) 


(14) 


The distributed voltage V(x,t) is the control applied to the system. Equation 
and the boundary conditions (15) form the control system model. 


IT 


3 SENSOR MODEL 

A distributed sensor is the one whose output is a function of structural responses 
at different locations. It can be a group of point sensors or a spatially continuous one. 
These responses are observed either discretely or continuously in space. Using the 
latter has the advantage that complicated computations based on point measurements 
can be reduced because the sensor geometry itself provides the processing. The spatial 
aliasing from an array of point sensors can be avoided. Typical noncausal sensor 
dynamics such as gain rolloff without phase shift is possible by using distributed 
sensors [3]. 

Figure 2 shows the sensor structure. PVDF is strain sensitive as it relies on the 
piezoelectric effect to produce the electric charge. The charge is proportional to the 
strain induced by the structure. This type of sensing is actually an inverse process of 
piezoelectric actuation. Based on the constitutive equation, the induced charge per 
unit length from the sensor strain is 


q(x,t) ~ ~Esd 3 x e s . 


(16) 
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PVDF layer 



Figure 2: The sensor model 


The sensor strain is related to the beam curvature by 

hf) + h s d 2 w 


e. = 


2 dx 2 ' 

The electrical charge along the beam is 

Q{x,t) - i q(x,t)F(x)dx 
J o 


= - E,d 


hb + h s [ x 


31 ' 


f F(x) 
Jo 


d 2 w 

Ih 2 


dx 


(17) 


(18) 


where F(x) is the weight function or shape function of the sensor. It is the local 
width of the electrodes covering both sides of the sensor layer. The function F(x) can 
be designed according to the need for interpreting the sensor signal. The capacitance 
between the electrodes of the sensor layer is 


CqC r As 

= h. 


(19) 


where £o and e T are the vacuum permittivity and relative permittivity constants 
respectively. The output voltage from the sensor is 


V.(x,t) 


Q{x,t) 

c 



( 20 ) 
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where 


( 21 ) 


js _ E s dzi(hb + h s ) 

s ~ 2 C 

is a constant. Suppose the sensor covers the whole beam, then 

fL d 2 u> 

V.(t) = -K.j Q F{x)—dx. 

Equation (22) is the sensor output equation. The output voltage is the weighted 
integration of the beam curvature. Integrating the output voltage by parts twice in 
spatial variable, we have another form of the sensor output, 


v,(t) 


dw(L,t) dF{L) 

-A, — )^F(L) + w(L,t)- 


dx 


dx 


K 


uj 


L i , 

w[x , t) — — - — d. x. 


dx 2 


>3) 


We shall see later from Equation (23) that different measurement outputs can be for- 
mulated to meet our control needs by choosing the appropriate sensor shape function 


4 DISTRIBUTED CONTROL ALGORITHM 

We design the control algorithm by Lyapunov’s direct method. The energy func- 
tion is used to measure the amount of vibration of the system. We need to find a 
control algorithm such that the closed loop system is asymptotically stable. One 
advantage of this method is that there is no need for model truncation. 

Given the system (14) with boundary conditions (15) and an energy functional 
E(t), we need to find a control V(x,t) such that 

lim E(t) = 0. (24) 

^—►oo 

It suffices to find a control V(x,t), such that 

<0, t > 0. (25) 

We define the energy function as follows: 



The first term is the stored energy due to bending. The second term is the kinetic 
energy due to the translation motion. The last term of the integrand is the kinetic 
energy from rotation of the beam element corresponding to Equation (14). a and b 
are positive constants. 
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Taking derivative of E(t) with respect to time and incorporating the system equa- 
tion (14) into it, we have 


dE{t) 

dt 


L 


L t d 2 w 

[a 


d 3 w 


+ 


dx 2 dx 2 dt 
c dw d 2 V 


d 2 w d 3 w I dw d 4 w 
+ b fadi'dxdt 2 + AlH'd^d? 
El dw d 4 w 


pA dt dx 2 pA dt dx 4 


}dx. 


Integrating by parts and setting 


a = 



We obtain 


dE(t) 

dt 


E7 cftxv dhv 
pA dx 2 dxdt 


r —L 


x — 0 


El d 3 w dw x ~ L 
pA dx 3 ' dt x=0 


I dw d 3 w 
A dt dxdt 2 


X ~ L c ft dw d 2 V , 


Introducing the boundary condition (15), we have 


dEjt) 

dt 


pA 


V(L,t) 


d 2 w 

dxdt 


I dw d 3 w 
A' dt dxdt 2 


C dV(L,t) dw 
pA dx dt ; 
c r L dw d 2 V 
+ pA Jo dt ' dx 2 


(28) 


(29) 


(30) 


The first term in Equation (30) contains the rotational velocity of the beam at the 
end. The second term has the force applied by the actuator. The third term is 
the product of the velocity of the displacement and the angular acceleration at the 
end of the beam. There is a second partial derivative of V(x,t) with respect to the 
spatial variable. We can design the appropriate modal controller by choosing the 
right V{x, £). Our purpose here is to find the control such that the time derivative of 
the energy function is negative. 

Let V(x,t ) be decomposed as the product of a spatial and a time function 


V{x,t) = u(s)g(0 


(31) 


where u(x) is the actuator shape function, q(t) is the coordinate function. We assume 
that the function u(x) has continuous second derivative on the interval (0, L) and has 
compact support over the interval; then the first two terms in Equation (30) vanish. 
Since the third term is negative from its physical meaning it will not cause energy 
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increase. We need to analyze the influence of the integration term in the equation. 
Substituting the voltage function into the last term in Equation (30) 


c r L dw d 2 V , c r L dw d 2 v(x) 

7 a L We^ dx= ^ {t) l m-e^ dx ( 

where q(t) is the time coordinate of the controller. 

We further introduce the feedback control by using the sensor output signal V 


32) 


q(t) 


dVsjt) 

dt 

r , f L dw(x,t)d 2 F(x) 

-'Wo St a^~ dz 


Then (32) becomes 

dEjt) 

dt 


< 


< 


c /4X f L dw d 2 v{x) 
jA q{t) Jo ' dt ' !h 2 


dx 


c f L dw d 2 F(x) f L dw d 2 v(x) 

* s pA Jo dt * dx 2 X Jo dt dx 2 


0 


(33) 


(34) 


Hence the system is asymptotically stable. The feedback control is given by 


V(x,t) 



dw(x, t) d 2 F(x) 
dt dx 2 


dx. 


(35) 


The introduced control is velocity feedback control. It takes into account the 
bending rate along the beam and introduces damping to the system. Here there is 
no need for the model modal truncation. 

When the control V(x ,t) is uniformly distributed in space, = 0, if we further 
assume that there is no elastic bonding layer to be present between the piezoelectric 
and the substructure, that is, there is no shear lag between the two layers, the strain 
is transferred between the piezoelectric and the beam over an infinitesimal distance 
near the end of the actuator [4]. We then have the simplified equation, 


EI dSv{x 1 t) + p/{ d 2 w(x,t) 


dx 4 


dt 2 


dSvj^t) _ 

pI dx 2 dt 2 


with boundary conditions 


u;(0, t) 
du;(0, t) 


dx 


El 

El 


d 2 w(L,t ) 

dx 2 

d 3 w(L, t) 

dx 3 


0 

0 


c.V(t) 


0. 


(36) 


(37) 
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V depends only on time in this case. This is a boundary control problem. We again 
use Equation (26) as the energy function. Repeating the procedure, we observe that 
only the first term in Equation (30) survives. This gives us 

L 

Notice that it is sufficient to choose 


to make 


dE(t) 

dt 


d 2 w 

(39) 

dxdt r 

x—L 

< 0 

(40) 


so as to asymptotically stable the system. The tip rotation speed is available from 
the sensor output Equation (23). 


5 SENSOR AND ACTUATOR SHAPE 

CONSIDERATION 

It is interesting to see that by introducing velocity feedback controller (35), the 
energy decay rate is given by (34) whose right hand side is a function of F{x) and 
v(x). Here F(x) and u(:r) are the shape or weight functions of the actuator and the 
sensor. They add weight for the control and measurement at each cross section along 
the beam. If we consider the displacement of the beam as the sum of a series of 
products of modal function and its coordinate, we can further analyze the effect of 
sensor and actuator shapes to different vibration modes. 

When the electric field is applied to the piezoelectric lamina, the actual piezoelec- 
tric actuation happens only in the region where both sides are covered by the elec- 
trodes. The same is true for collecting charge from the sensor layer. Hence, changing 
the width of the layout of the conductor is equivalent to varying the weighting func- 
tions. In this sense, it is possible to design the controller to suppress a particular 
mode or to design a distributed sensor to measure an interested mode. 

Consider the sensor output (23). We can get different information from the system 
by tailoring the right weighting function F(x). For example, we may select F(x) in 
such a way that 


d 2 F{x) 

dx 2 

dF(L) 

dx 

F(L) 


0, 0 < x < L 

0, x = L 

1 
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the integral term vanishes, so does the second term, and then the sensor output 


Vs(t) = 


c)w(Lj) 

dx 


■ 12 ) 


represents the angular deflection at the tip of the beam. Similarly, we can measure 
the displacement of the tip. From (20), the sensor can be used to build strain gauge 
by setting F(x) to be a spatial Dirac delta function 6(a;). 

We think that it is theoretically possible to use segmented sensors for the control 
of flexible structures. Digital control provides the ability to implement a sensing 
network with simple computation to rearrange the sensor layout and get different 
measurements with one sensor layer. Some measurements which are difficult to obtain 
in conventional way may be available by using distributed sensors. It may be feasible 
to implement full state feedback in relevant semigroup control formulations. 

The effectiveness of the PZT controller in introducing structural damping and its 
influence to the system dynamics of realistic size is based on the control authority ol 
the controller. Using velocity feedback shall increase the damping, but the control 
gain is limited to the electric field limit to avoid depolarization of the actuator. The 
actuator weighting function u(x) also plays a role here. We know that the bending 
moment is concentrated mostly at the end of the actuator of the beam. Hence more 
weighting should be placed on the region with high average strain. 

The feedback control (35) actually provides Voigt type damping since the rate of 
change of the bending curvature is used for feedback (22) and this rate is proportional 
to the rate of change of the structural strain. The augmented composite beam has an 
altered constitutive equation. The stress is no longer just proportional to the strain, 
but a linear combination of strain and the rate of strain change with respect to time. 

We now analyze the effect of both sensor and actuator shape functions to the 
damping control of different vibration modes. We use a Ritz-Galerkin procedure to 
implement modal expansion. We write the beam displacement w(x,t) as 


w{x,t) = ®k(x)dk{t) (-13) 

k = l 

where <J>^(:r) is the modal function and d k is the time coordinate. We can choose the 
orthogonal modal functions. We rewrite here the control form of the previous section 


V(x,t) = v{x)q(t). 


(44) 


Substituting the modal forms into Equation (14), multiplying each term with 4>j(.r) 
and then taking spatial integration along the beam, we get 


£ y $ t (x){pA$ k {x) - pl^\x))d k + EM l (x)^[ 4 \x)d k dxj = cq(t) J ^(x)iP\x)dx 

(45) 
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where d and d stand for the first and second time derivatives of the function d(t): 
stands for the ith spatial derivative of v(x). We then have 



rt n 

mikdk + kkdk = Qi 

k= 1 k - 1 

n 

(46) 

where 

Ql = C lkdk{t) 

k= 1 

(47) 


m i k = y] / ^/(x)[p/14>a(x) - p/$A (2) (x)] dx 
k= l J 

(48) 


kik = E / EI^(x)^\x)dx 

k=l J 

(49) 


Cik = cq(t) J ^i(x)v^ 2 ' l (x)dx. 

(50) 


The compact modal form is 

Md(t) + Kd(t) = Q(t) (51) 


where d(t) and Q(t) are nth order column vector functions. M is the inertial matrix, 
K is the stiffness matrix, Q(i ) is the modal control input. The damping of different 
elastic modes is influenced by 

Q(t) = -Cd[t) (52) 


and 

Md(t) + Cd{t) + I<d(t) = 0 (53) 

where C is a n by n damping coefficient matrix. Its elements are derived from the 
control law (35) 

d 2 F{x) 


Ctk 


= cK , [ 

Jo 


1 d 2 v(x) 

— r-br^-dx 


L 




-dx. 


(54) 


dx 2 Jo dx 2 

Observing Equation (54), we notice that in addition to the control authority de- 
termined by the actuation and sensing constants c and K a the added damping to 
a specific mode depends on the shape functions F(x) and v(x). The function i>(.r) 
in the first integral decides the amount of control effort applied to the /th mode. 
Similarly, F(x) provides the observation of the A:th elastic mode. The coefficient c/a- 
can be viewed as a measure of the damping to the /th mode by control based on the 
information from the kth elastic mode. If we choose the sensor shape to be such a 

function that 

d 2 F(x) 


dx 2 


$*(*), 


35) 


we can measure the A;th mode completely. When the second spatial derivative contains 
several modes, we shall get the combined information from the sensor. The similarity 
holds for the actuator, too. Properly selecting F(x) and v(x), we can observe and 
suppress the vibration modes. 
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6 CONCLUSIONS 


We have embedded a static PZT actuator model into the improved Euler-Bernoulli 
beam model to form a composite beam model with the rotational inertia effect con- 
sidered. We further used a distributed PVDF sensor to measure the elastic bending 
modes. A closed loop controller has been designed by using Lyapunov’s direct method. 
The closed loop system extracts energy from the system. The closed loop system is 
asymptotically stable. Finally, we have discussed the effects of different sensor and 
actuator shapes to the elastic modes. We point out that it is feasible to select suitable 
sensing and control weight to implement vibration control to some specified elastic 
modes. 

Further research is needed regarding aspects of estimation of the energy decay 
rate and real time implementation of the control law. We also would like to consider 
modeling the substructure with the Timoshenko model or the geometric exact rod 
model. Also, the real impact of the modal controller needs to be verified and further 
explored by experiments. 
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